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Abstract 

In this paper a class of discrete optimization problems with uncertain costs is discussed. 

The uncertainty is modeled by introducing a scenario set containing a finite number of 
cost scenarios. A probability distribution over the set of scenarios is available. In order to 
choose a solution the weighted OWA criterion (WOWA) is applied. This criterion allows 
decision makers to take into account both probabilities for scenarios and the degree of 
pessimism/ optimism. In this paper the complexity of the considered class of discrete 
optimization problems is described and some exact and approximation algorithms for 
solving it are proposed. Applications to a selection and the assignment problems, together 
with results of computational tests are shown. 

Keywords: robust optimization, weighted OWA, computational complexity, approximation 
algorithms 


1 Introduction 

Most practical decision making problems arise in a risky or uncertain environment, which 
means that an outcome of each decision is unknown and depends on a state of the world, 
which may occur with some positive probability. If probabilities for the states of the world 
are available, then each decision leads to a lottery, i.e. a probability distribution over the set 
of all possible outcomes. A decision problem can then be reduced to establishing an ordering 
of the set of lotteries. According to the classic expected utility theory by von Neumann and 
Morgenstern [zzmz!, the decision maker can assign an utility to each outcome, if he accepts 
some simple and appealing axioms. He can then compute an expected utility of each lottery 
and choose a decision which leads to a lottery with the largest expected utility. 

The expected utility can be seen as a weighted average of outcomes, where the weight 
of each outcome is just the probability of obtaining it. Thus, in the von Neumann and 
Morgenstern theory, the weights are independent of the outcomes and other probabilities 
of the lottery. However, it has been observed in human behavior that this assumption is 
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often violated (see [6] for a deeper discussion on this topic). Many decision makers pay more 
attention to unfavorable outcomes and would assign larger weights to such outcomes. In such 
a situation the weight of each outcome depends not only on its probability, but also on its rank 
in the lottery. Such weights may better reflect the pessimism/optimism of decision makers. 
A theory of such rank dependent, transformed probabilities was introduced by Quiggin [53] 
(see also [25]h 

In many practical situations the probabilities of scenarios are not available. We then obtain 
a decision problem under uncertainty. In this case, decision makers may assign subjective 
probabilities to scenarios m and compute the expected utility with respect to these subjective 
probabilities. However, determining the subjective probabilities may be not an easy task. An 
alternative approach is to apply some decision criteria such as the min-max, min-max regret, 
Hurwicz, or Laplace ones. In particular, in the Laplace criterion we apply the principle 
of insufficient reason and assign equal probability to each scenario. Each decision is then 
evaluated as the average utility of all possible outcomes. For a deeper discussion on decision 
making under uncertainty and description of the criteria we refer the reader to m- 

In this paper we discuss a class of discrete optimization problems, in which a finite set 
of feasible solutions is specified. In the deterministic case a cost of each solution is known 
and a decision problem consists in choosing a solution with the minimum cost. Discrete 
optimization problems are often represented as integer programming ones, in which the set 
of feasible solutions is described in compact form by a system of constraints. A class of 
deterministic discrete optimization problems was described, for example, in [ 21 ]. In many 
practical situations, the cost of each solution is unknown and depends on a state of the 
world which may occur with some positive probability. Each state of the world induces 
a cost scenario. A scenario set containing all possible cost scenarios is part of the input. 
In this paper we assume that this scenario set contains a finite number of explicitly listed 
scenarios. We also assume that probabilities for the scenarios are available. Notice, that 
under uncertainty, the principle of insufficient reason can be applied, which assigns equal 
probabilities to scenarios m- In order to choose a solution, we will apply the Weighted 
Ordered Weighted Averaging (WOWA for short) operator, proposed by Torra [ 26 ]. Given a 
solution, this operator allows us to define a rank-dependent weight for this solntion under each 
scenario. This weight can be seen as a distorted scenario probability and the WOWA criterion 
is then a special case of the Choquet integral with respect to distorted probabilities m- 
We can evaluate each solution as a weighted average of its costs over all scenarios. The 
WOWA criterion contains basic criteria used in decision making under risk and uncertainty, 
such as the expectation (weighted mean), maximum, minimum, Hurwicz, and Laplace ones. 
Furthermore, if the principle of insufficient reason is applied, then WOWA becomes the OWA 
criterion proposed by Yager [ 28 ] . 

If the uncertainty is represented by a discrete uncertainty set, it is common to use the 
robust approach m to compute a solution. In this approach we assume that decision makers 
are risk averse and we seek a solution which minimize the cost in the worst case. This leads to 
applying the min-max or min-max regret criteria to choose a solution. The traditional robust 
approach has, however, several drawbacks. The min-max criterion is extremely conservative 
and it is not difficult to show examples in which it gives unreasonable solutions m- In partic¬ 
ular, applying this criterion we may get a solution which is not Pareto optimal. Furthermore, 
the so-called drowning effect may also appear [7]. If the costs under some scenario are large in 
comparison with the costs under the remaining scenarios, then only this bad scenario is taken 
into account in the process of computing a solution (information connected with the remain- 
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ing scenarios is ignored). Hence, in many applications a criterion which takes into account 
all (or at least a subset) of scenarios is required. The traditional robust approach assumes 
also that no probabilities are available for the scenarios, which is not always true. By using 
the WOWA criterion we can overcome this drawback. We can use the information connected 
with scenario probabilities and soften the very conservative min-max criterion. Furthermore, 
the WOWA criterion is consistent with the theory of rank-dependent probabilities and, in 
consequence, can better reflect the real attitude of decision makers towards risk. This is 
particularly important when decisions are not repetitious, i.e. they are implemented only 
once. The WOWA operator allows us to establish a link between the stochastic and robust 
optimization frameworks. The distorted (rank-dependent) probabilities allows us to establish 
a trade-off between the expected and the maximum solution costs. 

In this paper we focus on the computational properties of the considered problem. Since 
the maximum criterion is a special case of the WOWA criterion, all negative results known 
for the robust min-max problems remain valid if the WOWA criterion is used. Unfortu¬ 
nately, the min-max versions of all basic discrete optimization problems become NP-hard 
even for two scenarios. This is the case for the shortest path, minimum spanning tree, min¬ 
imum assignment, minimum cut, or minimum selecting items problems mmm- All these 
aforementioned problems become strongly NP-hard and also hard to approximate when the 
number of scenarios is part of the input [laiiiiiii]. Furthermore, when the OWA operator 
is used to choose a solution, then network problems (the shortest path, minimum spanning 
tree, minimum assignment, minimum cut) are not at all approximable [15]. However, for an 
important case of nondecreasing weights in the OWA operator, there exists an approximation 
algorithm with some guaranteed worst case ratio and the aim of this paper is to generalize 
this algorithm to the more general WOWA criterion. In the existing literature, the OWA 
operator and the more general Choquet integral have been recently applied to some multiob¬ 
jective optimization problems in [T0l[9l[8|. In these papers the authors propose some exact 
methods for solving the problems, which are based on a MIP formulation and a branch and 
bound method. 

This paper is organized as follows. In Section O we present the problem formulation 
and show a motivation for using WOWA as a criterion for choosing a solution under risk 
and uncertainty. In Section [3l we recall some known complexity results for the considered 
problem. In Section we propose an approximation algorithm for solving the problem, 
which can be applied to a large class of discrete optimization problems. Section [5] describes a 
method of constructing a mixed integer programming formulation, which can be used to solve 
the considered problem exactly. This method will be adopted from m- Finally, in Section [G] 
we show applications of the proposed model to a selection and the assignment problems. This 
section also contains results of computational tests, which describe the efficiency of the MIP 
formulation and the quality of the solutions that are returned by the approximation algorithm 
designed in Section 01 

2 Problem formulation 

Let E = {ei, ...,€„} be a finite set of elements and let <I> C 2® be a set of feasible solu¬ 
tions. In a deterministic case, each element & E has a nonnegative cost q and we seek a 
feasible solution A G $, which minimizes the total cost F{X) = 'Yhe-ex^i- denote such 
a deterministic discrete optimization problem by V. This formulation encompasses a wide 
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class of problems (see, e.g., mm)- We obtain, for example, a class of network problems by 
identifying E with edges of a graph G and ‘h with some objects in G such as paths, spanning 
trees, matchings, or cuts. Usually, V is represented as an integer 0-1 programming problem 
whose constraints describe in compact form. 

Assume that the element costs are uncertain and their values depend on a state of the world 
which may occur with some positive probability. Each such a state of the world induces an 
element cost scenario (scenario for short) Cj = {cji ,..., Cjn)- Let scenario set L = {ci,..., Ck} 
contain K explicitly listed scenarios. Letp = (pi, ... ,Pk) be a vector of scenario probabilities, 
i.e. pj is the probability of the event that scenario Cj will occur. The cost of a solution X 
depends on scenario Cj E T and we will denote it by F{X,Cj) = Yletex^ji- Choosing a 
solution X leads to a lottery, i.e. a probability distribution {piF(X,Ci), ... ,pkF{X,Ck)) 
over the costs of X under scenarios in T. In order to choose the best solution, we need to 
evaluate each lottery. To do this, we should assign a weight to each scenario and compute 
the weighted average solution cost. Under the assumption that the weight of the jth scenario 
is equal to pj, we obtain that the weighted average is just the expected solution cost, i.e. 
f{X) = E[X] = Z,e[K]PjFiX,c,). 
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Figure 1: A sample Shortest Path problem with four scenarios Ci = (5,6,0,5,0), C 2 = 
(1,6,4,0,0), C 3 = (1,6,6,0,0), and C 4 = (2,6,6,0,0). The costs of all three paths under all 
scenarios are shown in the table. 

Consider the sample Shortest path problem depicted in Figure [TJ The set of elements 
E = {ei,..., 65 }, contains 5 arcs of network G and the set of feasible solutions consists of 
three paths Ai, A 2 , and A 3 from s to t in G. There are 4 costs scenarios with the probabilities 
0.5, 0.2, 0.2 and 0.1, respectively, and the costs of each path under these scenarios are shown 
in Figured! Path Ai = { 61 , 64 } has the smallest expected cost and thus should be chosen 
when the expected value is used. However, this choice may be unreasonable for some risk 
averse or pessimistic decision makers. Observe that the probability that the path Ai will have 
a large cost equal to 10 is equal to 0.5 which may be too large and cause some decision makers 
to reject Ai. On the other hand, the path A 3 = { 62 , 65 } has the smallest maximum cost and 
should be chosen when the min-max criterion is used and the probabilities of scenarios are 
ignored. Notice that the path A 3 has a deterministic cost equal to 6 . However, some decision 
makers may feel that path A 2 = { 61 , 63 , 65 } is better, since the probability that the cost of 
A 2 will be less than 6 is equal to 0.7 and the probability that A 2 will have a large cost, equal 
to 8 , is only 0.1. The sample problem illustrates that there is a need of defining aggregation 
weights, which would depend not only on the scenario probabilities, but also on the rank 
positions of the costs of a solution under scenarios. For example, risk averse decision makers 
would assign a weight larger that 0.5 to scenario Ci, when solution Ai is considered. 

Before we discuss such a criterion, which fulfills the above requirements, we recall an 
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aggregation criterion, called the Ordered Weighted Averaging operator (OWA for short) pro¬ 
posed by Yager in [28]. Let w = {wi,... ,wk) be a weight vector such that Wj G [0,1] for 
each j € [K], [K] ~ ^ denote the set Given a vector of 

reals a = (ai,..., ax), let fi be a sequence of [K] such that aa-(i) > • • • > ao-(i^)- Then 

owa(„)(a) = ^ 
j&[K] 


The choice of particular weight vectors w leads to well known criteria in decision making 
under uncertainty (see, e.g., m)- Indeed, if tci = 1 and Wj = 0 for j = 2,..., iL, then the 
OWA criterion becomes the maximum. If wk = 1 and wj = 0 for j = 1,..., K — 1, then it 
reduces to the minimum. More generally, if tCfc = 1 and Wj = 0 for j G [K] \ {k}, then the 
OWA criterion is the A:-th largest cost and, in particular, when k = [K/2\ + 1, then the A:-th 
largest cost is the median. If rcj = 1/K for all j G [AT], then OWA is the average. Finally, 
A wi = a and wk = 1 — a, for some fixed a G [0,1], and wj = 0 for the remaining weights, 
then the Hurwicz pessimism-optimism criterion is obtained. 

Using OWA it is not easy to take the probabilities of scenarios into account. In par¬ 
ticular, the expected value is not a special case of OWA. In the following, we will present 
an aggregation criterion, called the Weighted Ordered Weighted Averaging operator (WOWA 
for short) proposed by Torra |26j . This criterion generalizes OWA and allows us to define 
rank-dependent weights which also depend on scenario probabilities. 

Let V = {vi,... ,vk) be a weight vector such that Vj G [0,1] for each j G 
Eig[i^] Vj = ^- Let w* be a continuous nondecreasing function on [0,1], w* : [0,1] 

The domain interval [0,1] is partitioned by points 0 = ^<-^<-^<---< 

The function w* is linear on each subinterval 
function) and satisfies, for a given weight vector v, the 
w* ~ '^i<j ^ ^ Observe that w* is uniquely defined by v. In this paper we 

also make the assumption: vi > V 2 > ■ ■ ■ > vk- Thus w* is additionally a concave function. 
Figure [2] presents three sample functions w* for K = 5 for three weight vectors. 


[K] and 
-)■ [ 0 , 1 ]. 

j G [K], {w* is piecewise linear 
bllowing equations: u;*(0) = 0 and 



Figure 2: Three sample functions w* for K = 5 
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Let p = {pi,... ,pk) be an additional weight vector such that pj € [0,1] for each j G [K], 
'^j£[K]Pj ~ I- Given a vector of reals a = (ai,... ,aK), let ct be a sequence of [K] such that 
ao-(i) > • • • > Then the WOWA criterion is defined in the following way [26] : 

wowa(„_p)(o) = ^ 
iGl/c] 


where 

- W*C^Pa{i))- 

i<j i<j 

The value of Uj is a weight assigned to number Oo-q). It is not difficult to show that 
ojj G [0,1] for each j G [K], and ~ Figure [2] shows two boundary cases of the 

form of function w*, which are attained at = (1,0,... ,0) and = (l/A',... , l/LC), re¬ 
spectively. The vector models the weighted mean, i.e. in this case we get wowa („2 p)(a) = 
^j^[K]Pj^j- "Fhe vector models the weighted maximum, which in the case of uniform 
p = (l/IL,... , 1/iL) is the usual maximum operator. It is easily seen that for arbitrary 
V and uniform p = (1/iL,... , 1/iL), WOWA becomes the OWA operator. An easy com¬ 
putation shows that the WOWA operator is monotone, i.e. when a and a' are such that 
aj > a'j for all j G [K], then wowa(^p)(o) > wowa(„p)(a'). Since it is a convex combina¬ 
tion of the components of a, we have min^gj^] Uj < wowa^p(a) < max^^j^] aj. Additionally, 
J2j^[K]Pj^j — 'wowa(„p)(a) < maxjgjj^] aj holds, when ui > • • • > vk- 

We now apply the WOWA operator to the uncertain problem V and provide the inter¬ 
pretation of the vectors v and p. For a given solution X G let us define: 

WOWA(X) = wowa(„p)(F(X,Ci),.. .,F{X,Ck)). 

We thus obtain an aggregated value for X, by applying the WOWA criterion to the vector 
of the costs of X under scenarios in T. Given vectors v and p, we consider the following 
optimization problem: 

Min-Wowa V : min WOWA(X). 

The vector p = {pi,... ,pk) denotes just the probabilities for scenarios. The vector v 
models the level of risk aversion (or the degree of pessimism/optimism) of a decision maker. 
Namely, the more uniform is the weight distribution in v the less risk averse a decision maker 
is. In particular, = {1/K,... ,1/K) means that decision maker is risk indifferent and 
minimizes the expected solution cost. On the other hand, the vector = (1,0,... ,0) and 
the uniform vector p = (1/K,... ,1/K) mean that the decision maker is extremely risk averse 
and minimizes the solution cost assuming that the worst scenario for the computed solution 
will occur. In general, vector v allows us to model various attitudes of decision makers towards 
risk. Moreover, nonincreasing weights are consistent with the concept of robustness. Given 
a solution X, let a be such that F{X,c^(^y) > ■ ■ ■ > F{X,c^(^x))- Then, the value of iVj can 
be seen as a distorted, rank-dependent probability of scenario and WOWA(X) is the 

expected solution cost with respect to the distorted probabilities. Notice that Uj depends not 
only on the scenario probability but also on the solution X. 

Let us consider again the sample Shortest Path problem shown in Figured) Suppose 
that V = ( 0 . 5 , 0 . 3 , 0 . 2 , 0 ). The computation of the weights for paths Xi = {61,64} 

and X2 = {61,63,65} is shown in Figure [ 3 l For Xi we get F{Xi,ci) > F{Xi,Ci) > 
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F{Xi,C 2 ) > F{Xi,C 2 ,) and u = (0.8,0.08,0.12,0). Hence WOWA(Xi) = 0.8 • 10 + 0.08 • 
2 + 0.12 • 1 + 0 • 1 = 8.28. Observe that for Xi, the worst scenario Ci has the weight equal 
to 0.8, which is greater than pi = 0.5 and the best scenario C 3 has the weight equal to 
0, which is less than = 0.2. This example illustrates how the vector v distorts the 
scenario probabilities for solution Xi, by paying more attention to worse scenarios. In a 
similar way we compute the weights for path X 2 , obtaining uj = (0.2,0.36,0.44,0) and 
W 0 WA(X 2 ) = 0.2-8+0.36-7+0.44-5+0-5 = 6.32. Observe that W 0 WA(A: 2 ) < WOWA(Xi), 
so a risk averse decision maker would prefer solution X 2 over Xi, contrary to the case when 
the expected value is used as the criterion of choosing a solution. 



Figure 3: The weights wi,... ,^4 for paths a) Xi = { 64 , 64 } and b) X 2 = { 64 , 63 , 65 }. 


3 Complexity of the problem 

In this section we discuss the complexity of Min-Wowa V. Notice that Min-Wowa V 
becomes the Min-Owa V problem, discussed in m. when p = {1/K,... , l/K) and v is 
an arbitrary weight vector. If additionally v = (1,0,... , 0), then Min-Wowa V is the Min- 
Max V problem, widely discussed in the literature devoted to the robust discrete optimization. 
Hence all negative complexity and approximation results known for Min-Owa V and Min- 
Max V remain valid for MiN-WowA V. Let us recall that Min-Max V is usually NP-hard 
even when K = 2. In particular, this is the case for all basic network problems such as 
Shortest Path, Minimum Assignment, Minimum Spanning Tree, or Minimum Cut 
(see, e.g., Furthermore, when K is part of the input, then for all the aforementioned 

problems, Min-Max V is strongly NP-hard and also hard to approximate within any constant 
factor [H [ig. The problem complexity becomes worse when the maximum criterion is 
replaced with the more general OWA one. It has been shown in m, that all the basic 
network problems are then not at all approximable. This negative result holds when the 
vector V is arbitrary. However, for nonincreasing weights in v the following positive result is 
known: 

Theorem 1 f[T5]L When vi > V 2 > ■ ■ ■ > vk and V is polynomially solvable, then Min- 
Owa V is approximable within viK. 
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In the next section we will generalize Theorem [T] to Min-Wowa V. It has been shown 
in |15j . that Min-Owa V can be solved in pseudopolynomial time and even admits a fully 
polynomial-time approximation scheme (FPTAS), when K is constant and some additional 
assumptions for V are satisfied. We now show that the reasoning can be easily generalized 
to Min-Wowa V. Observe that wowa(„p)(a) is nondecreasing with respect to each aj in a. 
This fact immediately implies, that there exists an optimal solution X to Min-Wowa V, 
which is efficient (Pareto optimal), i.e. for which there is no solution Y such that F{Y,Cj) < 
F(X,Cj) for each j € [K] with at least one strict inequality. Notice also that each optimal 
solution to Min-Wowa V must be efficient when all components of p and v are positive. 
Thus it is sufficient to enumerate efficient solutions and pick up a one, say X*, with the 
minimum WOWA(X*). For some problems, for example when V is the Shortest Path or 
Minimum Spanning Tree problem, such enumeration of efficient solutions can be done in 
pseudopolynomial time, provided that K is constant, by using techniques given in [ 3 ]. Hence, 
for constant K, MiN-WowA Shortest Path and Min-Wowa Minimum Spanning Tree 
can be solved in pseudopolynomial time. 

In order to construct an FPTAS, we need a definition of an exact problem associated with 
V and scenario set F (see [H]). Given a vector {bi,... ,bK), we ask if there is a solution 
X € such that F{X,Cj) = bj for all j € [K\. Let us fix e > 0 and let Pe{^) be the set of 
solutions such that for all X G $, there is X € F(;(<I>) such that F{Y,Cj) < (1 -|- e)F{X,Cj) 
for all j € [K], Basing on the results obtained in [22], it was proven in [T8| that if the exact 
problem associated with V can be solved in pseudopolynomial time, then for any e > 0, the 
set F(;(<I>) can be determined in time polynomial in the input size and 1/e. This implies the 
following result (the reasoning is the same as in |15jL 

Theorem 2. If the exact problem associated with V can be solved in pseudopolynomial time, 
then Min-Wowa V admits an FPTAS. 

Proof. Let us fix e > 0 and let X be a solution with minimum value of WOWA(X) among 
all the solutions in P£(<k). From the results obtained in |18l 122] . it follows that we can 
find X in time polynomial in the input size and 1/e. Assume that X* is an optimal solution 
to Min-Wowa V. Define vector b* = ((1 -|-e)F(X*,ci),..., (1 + e)F{X* ,ck)). By the 
definition of P^i^), there exists a solution Y' € F’e(^) such that F{Y',Cj) < (1 -k e)F(X*,Cj) 
for all j G [K]. The choice of X and the monotonicity of WOWA implies WOWA(Y) < 
WOWA(Y') < wowa(„_p)(6*) = (1 -t- e)WOWA(X*). We have thus obtained an FPTAS for 
Min-Wowa V. □ 

It turns out that the exact problem associated with V can be solved in pseudopolynomial 
time for some particular problems V, provided that the number of scenarios K is constant. 
This is the case for Shortest Path, Minimum Spanning Tree and some other problems 
described, for example, in | 3 ]. However, it is worth pointing out that the running time of the 
obtained FPTAS’s is exponential in K, so their practical applicability is limited to very small 
values of K. In the next section we will construct an approximation algorithm, which can be 
applied for larger values of K. 

4 Approximation algorithm 

In this section we construct an approximation algorithm for MiN-WowA V under the as¬ 
sumptions that ui > U2 > • • • > vk and V is polynomially solvable. We will also assume that 


Pj > 0 for each j € [K]. When pj = 0 for some j G [K], then we can remove scenario Cj from 
r without affecting the problem. We first prove some properties of the WOWA operator. Let 
a = (oi,... ,aK) be a vector of nonnegative numbers. Let tt be any sequence of [K], Let us 
define 

/7r(o) ^ ^ 

i6[^] 

where Uj = w*P-K{i)) — ''J^*PTr{i)) w* is the piecewise linear function induced 
by the vector of weights v (as in the definition of the WOWA operator). Observe that 
/ 7 r(a) = wowa(„^p)(o) when the sequence vr is such that 0 ,^( 1 ) > • • • > a-K{K)- The following 
lemma expresses the intuitive fact that f-wia) is a lower bound on wowa(„p)(a). 

Lemma 1. Given any vector a = (ai,..., ax) and any sequence tt of[K]. Then wowa(^_p)(a) > 
fwia). 

Proof. Assume w.l.o.g. that ai > 02 > • • • > ax- Consider two neighbor elements 0 ,^( 2 ) and 
® 7 r(j+i) ™ ^ such that Let us interchange ajr(i) and a^(j+i) in vr and denote the 

resulting sequence by vr'. We will show that /jrK®) ^ firiO'), where the equality holds when 
®7r(i) = ®7r(i+i)- This will complete the proof since we can transform tt into a = ( 1 ,... ,K) 
by using a finite number of such element interchanges without decreasing the value of /jr and 
fa{a) = wowa(„^p)(a). It is easily seen that /^/(a) - U{a) = a;'a^(i+i) + - a;ia^(i) - 

uji+ia^(^i+i) = - Wi)a,r(i) - (wi+i - a;')a^(i+i). Equality to- + 0 ;'+^ = uji + Wj+i (see 

Figure 0^) holds, and so —uJi = oJi+i —uj[ = a. Hence fT^fa) —/jr(a) = Q;(ajr(i) — a 7 r(i+i))- 
Since w* is concave, we conclude that which yields a < 0 since 

p^(i+i) > 0. Hence /^/(a) > U{a). □ 



Figure 4: Illustrations of the proofs of Lemma [Hand Lemma [2j 


Lemma 2. Given any vector a = (oi,..., ax)- Then wowa(„p)(a) < viK Yljeix] Pj^j■ 

Proof. Since w* is concave and piecewise linear, it follows that for each j G 

[K] (seeFigurelKb)). Inconsequence, wowa(„_p)(a) = - '^j&[K]^^^P<Tij)aa{j) = 

Let Ci = wowa(„ p)(cij,..., cxi) be the aggregated cost of element ei G E over all scenarios. 
Let X be an optimal solution for the costs Ci, i G[n]. We begin with a general result: 
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Theorems. Given any X. T/ien WOWA(X) < iTt-i • WOWA(X). 

Proof. Let a be a sequence of [K] such that F(A,Cct(i)) > ••• > F{X and Uj = 
iYli<j P( 7 {i)) ~ iYli<j Pcr(i))■ The definition of the WOWA operator and Lemma [T] imply 

the following inequality: 

WOWA(X) '»«)• = E E S E 

3^[K] ei&X ei&Xj&{K] ei£X 

Using Lemma[2l we get q < viK 'f2j£[K]Pj^ji- Hence, from the definition of X, we obtain 

Ci < Y, Ci < Kvi Y, Y1 

ei£X Si&X ei£X j£[K\ 

Since ui > • • • > ux it follows that 

WOWA(A) > p,F(A,c,) =Y.P^Y. = E E P^^P' (3) 

j£[K] j^[K] ei£X ei£X j£[K] 

Combining ([T]), ([ 2 ]) and ([3]) completes the proof. □ 

Theorem [3] leads to the following corollary: 

Corollary 1. If vi > ■ ■ ■ > vk and V is polynomially solvable, then WOWA V is approx- 
imable within viK. 

The bound obtained in Corollary [1] is tight and the worst case instance for the approxi¬ 
mation algorithm is the same as the one shown in m- Observe that the approximation ratio 
depends on the weight distribution in v. The more uniform is the weight distribution the 
smaller is the approximation ratio. We get the largest approximation ratio equal to K, when 
WOWA is the weighted maximum. On the other hand, when vi = 1/K, i.e. when WOWA is 
the expected value, then we get an exact polynomial time algorithm for the problem. 

In many cases the deterministic problem V is NP-hard, but is approximable within a 
factor of 7 . In this case the following result can be established. 

Theorem 4. If vi > ■ ■ ■ > vk and V is approximable within 7, then Min-Wowa V is 
approximable within 'yviK. 

Proof. The proof is similar to the proof of Theorem [3l In order to get a solution for costs Cj 
a 7 -approximation algorithm is applied. It is then enough to modify inequality ([2]), so that 
EeiSX Ci < 7 YleidX Ci < iKvi Y.ei£X '^j£[K] Pj^P' The rest of the proof is the same. □ 

5 Mixed integer programming formulation 

In this section we design a mixed integer programming (MIP) formulation for Min-Wowa V. 
We will use the idea proposed in [20] (see also 01811191 for alternative formulations for the 
OWA operator). Let us associate a binary variable Xi G {0,1} with each element Ci G E. Let 
^ {0) 1 }" be the set of all characteristic vectors of <h. Each vector x = (xi,. .. ,Xn) G 
x(<I>) defines a feasible solution X such that Cj G A if and only if x* = 1. We will assume 
that x(<I>) can be described by a set of linear constraints involving variables xi,..., Xn. Prom 
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now on we will identify a feasible solution X G $ with the corresponding characteristic vector 
X € x(‘h). Let us fix a feasible solution x G Let a be such that -F(x,Co-(i)) > ■■■ > 

F{x,Cfj(^x))- Define vector a = (ag, ai,..., ax) such that Oj = Ylj<iPcr{j)^ ag = 0. Let 
us define hx{0) = F{x,Cfj(^i)) for ai-i < 0 < at, i € [K], 9 G (0,1]. The following equality 
holds [2U] : 

WOWA(x) = K Vj / ^ hxi9)de. (4) 

j€[K] 



Figure 5: Illustration of formula (jl]) for vector (8,5,7, 3). In (c) the right-hand edge of the 
triangle in (b) is scaled by the factor l/{Kv 2 )- 

Formula ([3]) is illustrated in Figure [5l Observe first that the value of WOWA(x) is equal 

to the size of the area of the grey rectangles in Figure [5^. In Figure [5)3 the portion of this 

area that touches the triangle with base [1/iL, 2/77] is shown. After scaling the right-hand 

edge of this triangle by the factor l/{v 2 K) we obtain the area shown in Figure [S]:. After 

rotating the rectangles we obtain the area shown in Figure [5]i. Now it is easy to see that the 

_ 2 _ 

size of this area is equal to /J/ hx{6)d6. Multiplying it by V 2 K, we get the size of the area 

n 

from Figure [ 5 ) 3 . 

Equality (|3|) has the following interpretation (see also |20]). The value of K hx{9)d9 

~K~ 

is the average within the jth portion of 1/77 largest solution costs. Then WOWA(s) can be 
seen as the value of the OWA operator applied to these averages. When pj = 1/77, j G [77], 

J_ 

then 77 hx{9)d6 = F{x,C(j(j)) and WOWA(s) becomes the OWA aggregation operator. 


II 



















































Let us rewrite ([H) as follows: 


WOWA(a;) = K 


vj ( h^{e)d9- [ "" /i*(0)d0^ = 




10 


/o 


K 

K 


K-1 


AT / hx{9)d6 - ^ Uj+i / hx{9)d0 


yj=i i=o 

Define vk+i = 0. Since hx{0)d6 = 0, we have 


K „± 

K 


K 

K 


WOWA(x) = K lY^j / hx{e)d9 - / hx{9)d9 


ii=i 


i=i 


and we get the following equality: 


WOWA(x) = K Y - ^i+i) / "" hxi9)d9. 

je[K] -^0 


(5) 


Let us denote Lj{x) = Jq^ hx{9)d9 and u'- = vj — Vj+i, j G [K], Observe that u'- > 0 for all 
j € [AT + 1], by the assumption that vi > V 2 > ■ ■ ■ > vk- We are now ready to design a MIP 
formulation. In order to do this we adopt the idea from [20]. Observe first that the value of 
Lj{x) for a fixed x can be computed by solving the following linear programming problem: 

max 'Eke[K]^kF{x,Ck) 

'^kG[K] ^k = (6) 

d<Zk<Pk k£[K] 

Indeed, Ljix) can be computed in a greedy way. Let a be such that > ••• > 

F{x,c^(k))- We first allocate to the interval [0,j'/iL] the largest possible portion of Pa(i), 
then the largest possible portion of Pa{ 2 ) until [0,j'/iL] is completely filled. This is 

equivalent to solving (|6|). The dual to (|6|) for a fixed x and j takes the following form: 

min iPj+ J2ie[K]Pi(^ij 

f3j + aij > F{x,Ci) i G [K] (7) 

otij >0 z G [K] 

The strong duality theorem implies that Lj{x) equals the optimal objective value of ([7|l. 
Using ([5|) and d?]) we get that Min-Wowa V is equivalent to the following problem: 

min K ■ + E*e[A] 

13j + aij > F(x,Ci) ie[K],j e [K] 

aij >0 z G [K],j G [K] 

{xi, ...,Xn)€ x(^) 

We obtain a MIP formulation by substituting F{x,Ci) = Y2k€[n]XiCik and replacing the ex¬ 
pression (xi,..., Xn) G with a system of linear constraints involving xi,..., Xn- In the 

next section we will apply the MIP formulation to a sample problem. 
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6 Computational tests 


In this section we present the results of some computational tests. The tests were performed 
for two particular discrete optimization problems, namely the Selection and Assignment 
problems. We wish to verify the following two questions; 

1. How efficient is the MIP formulation, i.e. how the computation time, required to solve 
the MIP model, depends on the number of elements n in the set E, the number of 
scenarios K and the weight distribution in v? 

2. What is the quality of the approximation algorithm designed in Section SP 

For both problems we used the following method of generating the tested instances. For 
each scenario j E [K] we chose aj, which is a random integer uniformly distributed in [1,100], 
and then set pj = Oj/J ^ 1^]^ obtaining a positive probability for each scenario. 
In order to fix the weights vi,...,vk for scenarios, we used generating function ga{z) = 
Y^(l — a^) where a E (0,1) is a fixed parameter. Notice that ga(z) is concave and is such 
that ga(0) = 0, ga(l) = 1- Given a and K, we set Vj = ga{j/K) — ga{{j — ^)/K) for j E [K] 
(see Figure [6]). The value of a expresses an attitude of the decision maker towards risk. The 
smaller the value of a, the more risk averse the decision maker is (the less uniform is the 
weight distribution vi,... ,vk)- For each generated instance the CPLEX 12.5 solver with 
standard settings was used to solve the corresponding MIP formulation. We fixed the time 
limit to 3600 seconds. The solver was executed on a computer equipped with a 2.5 GHz 
processor with 8 GB RAM. 



Figure 6: Generating function ga{z) for a E {10 ^,10 ^,10 ^,10 and the weights 
ui,..., U 4 for A = 4 and a = 10“^. 


6.1 The selection problem 

In this section we apply the MIP formulation and the approximation algorithm designed in 
Section 0] to the following Selection problem. Assume that A is a set of n items and we wish 
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to choose exactly q of them to minimize the total cost. Hence $ = {X C E : \X\ = q}. The set 
of characteristic vectors x(‘h) can be described by one constraint of the form xi + - ■ -Ezn = q, 
where xi, ... ,Xn are binary variables associated with the items in E. The Selection problem 
has been recently discussed in a number of papers. Its min-max version has been proven to be 
NP-hard for two scenarios [3] , strongly NP-hard and hard to approximate within any constant 
factor when the number of scenarios is part of the input m- Hence the same negative results 
hold for Min-Wowa V. 

We performed the tests for the number of items n chosen from the set {160,200}, the 
number of scenarios K chosen from the set {2,... , 20}, and the parameter a chosen from the 
set {10“^, 10“^, 10“^}. We also fixed q = 0.25n, i.e. we assumed that exactly 25% of the 
items must be chosen. Under each scenario the cost of item is an integer that is chosen 
randomly with uniform distribution from the set {0,..., 100}. For each combination of n, 
K and a we have generated 10 random instances. We first applied the MIP formulation to 
obtain the optimal solutions for the instances. The computational times required by CPLEX 
to solve them are shown in Figures [3 and [HI It can be observed that the computational times 
quickly grow with the number of scenarios. The problem is also harder to solve for smaller 
values of a. We were unable to solve any instance with n = 200, K = 20, and a = lO”"^ 
within the time limit of 3600 s. 



a= 10-2 



Figure 7: Computational results for Min-WOWA Selection - running times for n = 160 
and all combinations of a and K. The solid line shows the average computational time. The 
numbers in brackets show the number of instances which were not solved within 3600 s. 

We next applied the approximation algorithm, designed in Section jH to the generated 
instances. The obtained results are shown in Figures P and [TOl 

For each instance, for which an optimal solution was known, we computed the percentage 
deviation of the cost of the approximate solution from the optimum. The quality of the 
solutions returned by the approximation algorithm seems to be good in comparison with the 
worst theoretical performance. The largest reported deviations from optimum are not greater 
than 32%, whereas the largest theoretical deviation varies from about 90% (for K = 2 and 
a = 10“^) to about 638% (for AT = 20 and a = 10“^). We can thus obtain reasonable 
solutions as long as the distribution of the costs under scenarios is uniform. It is interesting 
that the deviation from the optimum depends more on a than on K. The performance of 
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Figure 8: Computational results for MiN-WOWA Selection - running times for n = 200 
and all combinations of a and K. The solid line shows the average computational time. The 
numbers in brackets show the number of instances which were not solved within 3600 s. 


the approximation algorithm is clearly better for larger a, i.e when WOWA is closer to the 
expected value. On the other hand, for a fixed a, the performance is significantly better only 
for K = 2. For AT = 4 it becomes worse. Interestingly, one can observe a slightly better 
performance when K increases from 4 to 20. 

6.2 The assignment problem 

In this section we apply the MIP formulation and the approximation algorithm designed in 
Section |4] to the following Assignment problem. We are given a bipartite network G = 
(Fi U V 2 ,E), where n Pa = 0, Itil = 11^21 = m, and E = Vi x V 2 , jEj = n = ml Set 4> 
contains all subsets of E which form a perfect matching (assignment) in G. We can associate 
a binary variable Xij with each element (edge) Cij € E and the set of characteristic vectors 
x(^) can be described by the assignment constraints of the form X]iG[m] ~ ^ J ^ I™"] 
and Ylje[m] ^ij ~ ^ ^ [^]- min-max version of the Assignment problem is known 

to be strongly NP-hard and hard to approximate within any constant factor m- 

We performed the tests for the number of nodes m (|Pi| = IV 2 I = m) chosen from the set 
{40,50}, the number of scenarios K chosen from the set {2,... ,20}, and the parameter a 
chosen from the set {I0“^, 10“^, 10““^}. Observe that the cardinality of E was 1600 and 2500, 
respectively. Under each scenario the cost of element is an integer that is chosen randomly 
with uniform distribution from the set {0,..., 100}. For each combination of m, K and a 
we have generated 10 random instances. We first applied the MIP formulation to obtain the 
optimal solutions for the instances. The computational times required by CPLEX to solve 
them are shown in Figures [11] and [T2j One can notice that the computational times quickly 
grow with the number of scenarios. The problem is also harder to solve for smaller values 
of a. For AT = 10 and m = 50 some instances could not be solved within the time limit of 
3600 s. 

We next applied the approximation algorithm, constructed in Section H] to the generated 
instances. The computational results are presented in Figures [T3| and dH We can derive 
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Figure 9: The performance of the approximation algorithm run on the instances of Min- 
WOWA Selection ~ percentage deviations from the optimum for n = 160 and all combi¬ 
nations of a and K. The solid line shows the average deviation. 


32 . %Err 
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a = 10-2 



Figure 10: The performance of the approximation algorithm run on the instances of MiN- 
WOWA Selection - percentage deviations from the optimum for n = 200 and all combi¬ 
nations of a and K. The solid line shows the average deviation. 


similar conclusions as for the Selection problem. The largest deviation from optimum 
reported was about 30%, which is much less than the worst theoretical performance which 
varies from about 90% (for K = 2 and a = 10“^) to about 501% (for A = 10 and a = 10“^). 
As for the Selection problem, the deviations from optimum is significantly smaller only for 
K = 2. 

7 Conclusions 

In this paper we have discussed a wide class of discrete optimization problems in which the 
uncertain costs are specified in the form of a discrete scenario set. A probability distribution 
over this set of scenarios set is provided. We have applied the weighted OWA criterion to 
choose a solution. This criterion allows us to take both scenario probabilities and attitude of 
decision makers towards a risk into account, as the weights assigned to scenarios are distorted 
(rank dependent) probabilities. Our approach contains the traditional robust (min-max) and 
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Figure 11: Computational results for MiN-WOWA Assignment - running times for m = 40 
and all combinations of a and K. The solid line shows the average computational time. 


stochastic approaches as special cases. The problem of minimizing the WOWA criterion is typ¬ 
ically NP-hard for two scenarios. It becomes strongly NP-hard and also hard to approximate 
when the number of scenarios is part of the input. It is thus important to provide efficient 
approximation algorithms for the problem. One such an algorithm has been constructed in 
this paper. It can be applied, if the underlying deterministic problem is polynomially solv¬ 
able. The efficiency of the MIP formulation and the quality of the approximation algorithm 
were tested for two particular problems, namely the selection and the assignment problems. 
The MIP formulation can be used when the number of scenarios is small. For larger number 
of scenarios the approximation algorithm may be an attractive choice. The performance of 
the approximation algorithm seems to be good when the element costs under scenarios are 
chosen randomly. It may be poorer for more correlated costs and investigating the quality of 
the algorithm in this case requires additional tests. 
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